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A DYNAMIC VISCOELASTIC ANALOGY 
FOR FLUID-FILLED ELASTIC TUBES 

ANDREA GIUSTL AND FRANCESCO MAINARDI^ 


Abstract. In this paper we evaluate the dynamic effects of the fluid viscosity for fluid 
filled elastic tubes in the framework of a linear uni-axial theory. Because of the linear 
approximation, the effects on the fluid inside the elastic tube are taken into account 
according to the Womersley theory for a pulsatile flow in a rigid tube. The evolution 
equations for the response variables are derived by means of the Laplace transform tech¬ 
nique and they all turn out to be the very same integro-differential equation of the 
convolution type. This equation has the same structure as the one describing uni-axial 
waves in linear viscoelastic solids characterized by a relaxation modulus or by a creep 
compliance. In our case, the analogy is connected with a peculiar viscoelastic solid which 
exhibits creep properties similar to those of a fractional Maxwell model (of order 1/2) for 
short times, and of a standard Maxwell model for long times. The present analysis could 
And applications in biophysics concerning the propagation of pressure waves within large 
arteries. 


1. Introduction 

It is well known that the dynamic theory for fluid filled elastic tubes might find ap¬ 
plications in the analysis of the propagation of pressure waves within large arteries. In 
the framework of a uni-axial theory, the effect of the viscosity of the fluid inside the tube 
[e.g. blood) is usually taken into account including a friction term in the one-dimensional 
momentum equation, due to the Poiseuille’s formula for steady flows in rigid tubes, see 
e.g. [12], and [2], 

Restricting ourselves to a linear approach, here we hnd it more appropriate to evaluate the 
friction term for a non-stationary flow, which is easily derived starting from a sinusoidal 
flow, as in the renowned Womersley model for pulsatile flow of viscous (Newtonian) fluids 
in rigid tubes [21], see also [5], [11]. 

In Section 2 we follow the approach displayed in [3] considering the Navier-Stokes equa¬ 
tions for an incompressible Newtonian fluid in cylindrical coordinates. For such system of 
PDEs we neglect the motion in the circumferential direction and we assume the flow to be 
quasi-one-dimensional. In particular, the radial velocity is assumed to be very small with 
respect to the axial one. The resulting equations are then integrated over the radial coor¬ 
dinate, reducing the number of independent variables to two (i.e. time t and axial distance 
x) and the dependent variables to three (pressure, averaged velocity, and cross-sectional 
area). A pressure-area relation for a uniform elastic tube, longitudinally tethered, is then 
used to eliminate one of the remaining dependent variables. However, the hnal equations 
contain two quantities which depend on the velocity prohle. 

In Section 3, working within the framework of a linear uni-axial theory, we reduce the 
system of PDEs discussed in Section 2 to a system of two linear equations that contains 
only an additional unknown function, i.e. the friction term, that we denote by fo{x,t) 
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and which still depends on the velocity profile. Then, because of our linearity assumption, 
we find it reasonable to evaluate /o in terms of the averaged velocity. This can be done by 
assuming the velocity profile to be described as in the case of an unsteady flow in a rigid 
tube with an unperturbed radius based on the Womersley theory for pulsatile flow. As a 
consequence, we are able to obtain the friction term corresponding to any time history of 
the pressure-gradient by means of a time convolution integral between a certain memory 
function <h(t) and the pressure-gradient itself. Finally, we obtain the evolution equation 
for the pressure waves that appears to have the same structure as the one known for stress 
waves in a viscoelastic solid (in the relaxation representation), see e.g. [14]. 

In Section 4 we discuss in detail this dynamic analogy in the framework of a linear 
theory of viscoelasticity. For further details about mathematical and historical aspects 
of this theory we recommend the interested reader to refer to specialized articles and 
treatises, e.g. [18], [13], [14], [19] and [8]. 

The memory function 4>(t) is then shown to be a complete monotonic function expressed 
in terms of a Dirichlet series. Therefore, the latter has to be intended as the rate of the 
relaxation modulus of a peculiar viscoelastic model. We then calculate the corresponding 
rate of creep T (t), which enables us to obtain the wave equation in the creep representation. 
Finally, by a suitable time integration of 4>(t) and 4'(t), we derive the corresponding 
relaxation modulus G{t) and the creep compliance J(t), providing a full characterization 
of the viscoelastic analogy. 

In Section 5 we derive the asymptotic representations of the memory functions 4>(t) 
and 'l'(t) for short and long times. The resulting functions show that the corresponding 
viscoelastic model is governed by a stress-strain relation of Maxwell type with a fractional 
order changing from 1/2 to 1 as time evolves from zero to infinity. For sake of clarity, we 
show the plots versus time of the functions 4>(t), 4'(t) and of G{t) and J{t). 

Finally, in Section 6 we complete the paper with concluding remarks and hints for future 
research. Furthermore, a detailed mathematical description of the technical results can be 
found in the Appendix. 


2. The basic non-linear equations 

Let us consider in a cylindrical coordinate system the Navier-Stokes equations of motion 
for an incompressible Newtonian fluid. If the motion in the circumferential direction is 
neglected, these equations read (denoting the partial derivatives with subscripts) 

( Px f Ur \ 

Ut + WUr + UUx H- = V [Urr ^ -h Uxx , 

12 I) \ p \ r J 

^ I Pr f Wr W \ 

Wt + WWr + UWx H-= U Wrr H-^ + Wxx , 

f p \ r J 

where t is the time, x, r denote the axial and radial directions, u, w the corresponding 
components of the fluid velocity and p the pressure. The constants p and u are the density 
and kinematic viscosity of the fluid. 

Then, following the approach shown in [3], the flow is assumed to be quasi one dimen¬ 
sional, i.e. the radial velocity w is very small with respect to the axial velocity u and 
the resulting equations are integrated over the radial coordinate r. As a consequence we 
arrive at Eqs. (14) and (15) in [3], which in our notation read: 
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( 2 . 2 ) 


Ut + ^(1 — x)^i + XUUx + — — 

A p K 

A, + 1(C/A)^0, 


R 




[ 


R 

2ru dr 


In these equations 

(2.3) 

is the averaged (uni-axial) velocity (the quantity measured in most experiments), 

(2.4) A = ttR^ 
is the cross-sectional area of the tube, and 

1 




2ru^ dr. 


^ t/2 Jo 

As a consequence, {U{x,t), A{x,t)}, the basic dependent variables for a quasi-one dimen¬ 
sional flow, tnrn out to be related also to the pressure gradient A = px, to the parameter 
X and to the friction term defined by 


( 2 . 6 ) 


= ^Ur 


R 


As a matter of fact, both quantities depend on the variation of the axial velocity u with 
r (the velocity profile) so they could not be properly determined without carrying ont a 
full two-dimensional non-linear treatment of the problem. It is worth stressing that when 
the fluid is inviscid the velocity profile is flat (x = 1) and the friction term is absent 
if{x,t) = 0). 

According to [3], in the presence of viscosity, a particnlarly simple assumption is the 
parabolic form of the velocity profile: 

(2.7) u = 2U{1 - r^/R^) , 
so that a = 4/3 and 

(2.8) f{x,t) =-8vU/R^ . 

This assumption, however, is too simple because it is satisfied only for stationary flows 
related to the Poiseuille approximation. 

3. The basic linear equations 
Linearizing Eqs. (2.2) we get the following systems of PDEs: 

Ut +Px/p = fo{x,t ), 


(3.1) 

from which we get 

(3.2) 


+ ^oUx — 0 , 


2^oUt + CqAx — Aq fo{x, t ), 

At + Aq Ux = 0, 


where Aq = irRn is the unperturbed cross-sectional area of radius Rq, 


(3.3) 


Cn = 


Ao 


dp 

Ja 


0 







4 


ANDREA GIUSTli AND FRANCESCO MAINARDI^ 


denotes the corresponding Moens-Korteweg velocity (for pressure waves in elastic tubes 
filled with an inviscid fluid) and 


(3.4) 


f 2Z/ 

Jo = -^Ur 
no 


Ro 


is the corresponding friction term due to the viscosity, evaluated at the same unperturbed 
radius. 

In the case of a parabolic velocity profile (Poiseuille approximation of steady flow), 
linearizing Eqs. (2.7) and (2.8), we get 

u = 2U{x,t) (1 -rVi?o) 


(3.5) 


fo{x,t) = - 


8U{x, t) 


T = 


Rl 


We note that the time constant r, introduced here for steady flow, would provide a proper 
time scale also in all profiles occurring for any unsteady flow. As a consequence, in every 
function of time appearing in our approach the constant r is to be interpreted as a proper 
time scale introduced by the viscosity of the fluid; in future, even when we will omit r in 
the time dependence, this scale-parameter is understood to be equal to one. 

In order to eliminate the velocity component in the radial direction, in the framework 
of a linearized theory for unsteady flow, and evaluate the corresponding friction term fo 
in terms of the uni-axial velocity U{x,t), we find it reasonable to take profit of the well 
established Womersley linear model for pulsatile flow of a viscous (Newtonian) fluid in 
rigid tubes [21]. 

Because of the linearity, we can easily compute the friction term corresponding to a 
pressure-gradient which is sinusoidal in time and oscillates with a frequency to, 

X{x,t]uj) = Ao(x) exp(—zwt) 

as derived by the Womersley model. Introducing the Womersley parameter^ 


(3.6) 


a := Ro\l — = 


/(XT , 


we get the friction expressed in terms of {^/uJt): 

(3.7) fo{xR]x) = S(a;) Ao(x) /p, 

with 

2 Ji{i^ 


(3.8) 


4'(a;) = 


where Jq, Ji denote the Bessel functions of order 0, 1, respectively, see e.g. [Ij. 

Then, being interested to transient motion, it is worth to pass from Fourier transform 
to Laplace transform with parameter s = ix so that we get 

2 /i(V^) 


(3.9) 


$(s) = 


^ /o(V^) ’ 

where /q, R denote the modified Bessel functions of order 0, 1, respectively. We have 
used the known relation between Bessel and modified Bessel functions of any order i/, 
JRiz) = Rluiz), see e.g. [Ij. 


^In our analysis we have preferred to express the Womersley parameter a in terms of the time scale t 
introduced in Eq. (3.5). 
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As a consequence we can formally obtain the friction term corresponding to any time 
history of the pressure-gradient A(x, t) by means of a Laplace convolution integral. Indeed 
we get: 

(3.10) fo{x, t) = d>(t) * A(x, t)/p 
from which we derive 

(3.11) Utix, t) = -[1 - $(f) *] A{x, t)/p 

where <h(t) denotes the inverse Laplace Transform of ‘h(s) and * denotes the time Laplace 
convolution^. 


So doing we have properly modified the relation between acceleration Ut and pressure 
gradient A for inviscid flow (the Euler equation) in terms of a Laplace convolution with 
a characteristic function <I>(t) in order to take into account the memory effects of the 
viscosity. We agree to call ‘h(t) the relaxation memory function, whose meaning will be 
clarified in the following with respect to our wave process. 

Indeed, denoting with Y = Y{x,t) a generic response variable (such as {U, A, p,A}), 
after simple manipulations we get the evolution equation 

(3.12) Yu{x, t) = Cq [1 - T>(t) *] ya,a,(x, t) 

which, in the absence of viscosity {foix,t) = 0, i.e. <I>(t) = 0), reduces to the classical 
D’Alembert wave equation 

(3.13) Ytt{x, t) = cl Y^^{x, t ). 

As a matter of fact, Eq. (3.12) is an integro-partial differential equation, which could 
be taken into account in order to investigate wave processes in fluid filled compliant tubes 
of biophysical interest. In particular, this evolution equation could be solved for x ^ 0, 
t ^ 0 together with a known input condition Y{0,t) = Lo(t) at x = 0, assuming that 
the tube is initially quiescent, and that there are no waves coming from x = -l-oo. This 
initial boundary-value problem is usually referred to as a signalling problem. The above 
simplifying assumptions are expected to explain the qualitative features of the distortion 
(due to the viscosity of the fluid) of the transient waves generated at the accessible end of 
the tube as they propagate along the tube. 


4. The dynamic viscoelastic analogy 


Erom Eq. (3.12) we easily recognize an analogy with the wave equation of linear vis¬ 
coelasticity in the Relaxation Representation, see [14]. For this purpose we will prove that 
the memory function <h(t) is completely monotonic (i.e. it is a non-negative, non-increasing 
function with infinitely many derivatives alternating in sign for t ^ 0), as required to rep¬ 
resent (with opposite sign) the non-dimensional rate of the relaxation modulus G{t) scaled 
with its initial value 0 < G(0’^) < oo according to 


(4.1) 


m 


1 dG 
G{0+) H ’ 


^We recall the definition of the time convolution between two locally integrable functions f(t), g(t): 
f(t)*g{t) := f f{T)g{t- T)dT = f f{t-T)g{T)dT, 

Jo Jo 

so that its Laplace transform reads f{s)g{s). 
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and consequently in Laplace domain 


(4.2) 


sGis) = G(0+) 1 - 5(s) 


We recall that G{t) represents the stress response to a unit step of strain, see [14] and it 
is usually assumed to be completely monotone as well. Being G(0'^) finite and positive, 
the corresponding viscoelastic solid exhibits stress waves with a finite wave front velocity 
CO = VG(0+)/p. 

In fact, carrying out the inversion of the Laplace transform (3.9), by means of the 
Bromwich formula and applying the residues theorem, we get the memory function in 
terms of a convergent Dirichlet series as 

4 °° 

(4.3) $(t) = - 2 

^ n=l 


where Xn are the zeros of the oscillating Bessel function Jq on the positive real axis. 
Further details about the proof can be found in the Appendix. So, we recognize that 
the relaxation memory function, resulting from the sum of a convergent Dirichlet series 
with positive coefficients, is indeed a complete monotone function, ensuring the complete 
monotonicity of the corresponding relaxation modulus G(t), with G(+oo) = 0. 

In Linear Viscoelasticity, the Relaxation Modulus G{t) is coupled with another material 
function, the Creep Compliance J{t), which represents the strain response to a unit step of 
stress. Hence, to complete the dynamic viscoelastic analogy we now consider, in addition 
to the relaxation representation of the wave equation (3.12), also the corresponding creep 
representation. 

We first observe that the two material functions G{t) and J{t) are connected through 
their Laplace transforms by the following fundamental equation 


(4.4) 


s J{s) 


1 

s G{s) 


that implies in the time domain 


(4.5) 


J(0+) = 


G(0+) 


^ 0 , J (-j-cx)) — 


G(+oo; 


■ = +00 


This means, according with the classification of viscoelastic media (see [14]), the model of 
our concern can be understood as of Type 11. 

Introducing the non-dimensional rate of creep 4'(t), defined as 


(4.6) 

we have 

(4.7) 


T(t) : = 


1 dJ 
J(0+) lit ’ 


sJ{s) = J(0+) 1 + T(s) 


Then, in the Laplace domain, the functions <I>(t) and 'l'(f) turn out to be related as 

r ~ 1-1 


(4.8) 


1 + 4'(s) 


1 - 4>(s) 


as a consequence of Eqs. (4.2), (4.4) and (4.7). 

Now, recalling the recurrence relations for the modified Bessel functions of the first 
kind, Iu{z), i.e. 
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with u = 1, after a simple manipulation we get: 


(4.9) 


$(s) 


2 hiV^) 


The inversion of the Laplace transform can be carried out using again the Bromwich 
formula. Then, applying the residues theorem we get: 


(4.10) 


S A ^ 

^(^) = - + - 2 exp (-/i^t/r) 
n=l 


where /in are the zeros of the oscillating Bessel function J 2 on the positive real axis. 

Therefore, we recognize that the 'll function, resulting from the sum of a convergent 
Dirichlet series with positive coefficients, is also a complete monotone function like 4>(t). 
As a consequence, the creep compliance J{t) turns out to be a Bernstein function (that is 
a non negative function with a completely monotone first derivative). We note, however, 
the presence of a positive constant term added to the Dirichlet series. This implies that 
the asymptotic expression for the creep compliance J(t), for long times, contains an extra 
term which is linear in time. The latter is a typical feature of the classical Maxwell model 
of viscoelasticity. 

Using the analogy with the linear dynamical theory of viscoelasticity, the wave equation 
corresponding to the creep representation reads, see [14], 

(4.11) [i + ^(t)*]ytt = cgy,,. 


where T(t) can be referred to as the creep memory function. Given the previous results, 
we can now trivially compute the Relaxation Modulus G{t) and the Creep Compliance 
Jit). 

In view of Eqs (4.1) and (4.6), we get G(t) and J(t) 


(4.12) 


Git) = G(0^ 


1 - 


f 


4>(t') dt! 


(4.13) Jit) = J(0+) 
Integrating term by term we get 

(4.14) 


1 + 


f T(t') dt! 

Jo 


n=l •- 


(4.15) 


G(t) = G(0+)<jl-2 

I 


A _ A e-A^t/r 


J(t) - J(0+) |l + — + ^ 


n=l ^ 


2 2 ® 




Due to the results for the the zeros of Bessel functions shown in [20], i.e. 

^ = 1, Jo(A„) = 0, Vne N, 


(4.16) 

00 

E 

n^l 

(4.17) 

00 

E 

n=l 


^2(/^n) = 0, Vne N 
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we get 


(4.18) 


Git) = G(0^ 


' 2 T2 (-'^nV'r) , 

n=l 


(4 81 °° 4 

(4.19) J{t) = + -2 (-/^n^/'T 

( ^ n=l 

As we see from Eq.(4.14), together with the condition G(+oo) = 0, we recover the 
identity in Eq. (4.16). This result gave us a relevant hint to derive Sneddon’s result about 
series involving zeros of Bessel functions of the hrst kind, by means of a technique based 
on the Laplace transform (see [7]). 



5. Asymptotic and Numerical results 

In this section we exhibit the plots (versus time) of the relaxation and creep memory 
functions as computed taking a suitable number of terms in the corresponding Dirichlet 
series (4.3) and (4.10), respectively. Erom our numerical experiments, we found that 100 
terms are surely suitable to get stable and confident results. 

We also compare these plots with the asymptotic representations of the two functions 
for short and long time in order to check their matching with the numerical solutions. Ac¬ 
cording to the Tauberian theorems, the asymptotic representations of <1>(1) and 'k(t), for 
short times, are formally derived by inverting their Laplace transforms (3.7) and (4.9), ap¬ 
proximated as ST ^ 00 respectively, see Appendix. Eor large times, the asymptotic repre¬ 
sentations are promptly obtained by taking the hrst term of the corresponding Dirichlet 
series. We get: 


(5.1) 


with Af 5.78 and 


m 



- exp (—Aft/r), t 


0 , 

00 , 


(5.2) 


4/(t) 



8 4 / 2 / \ 

- -h - exp (-//iVr), 
V r T 


t 


t 


0 , 


00 , 


with /if 26.37. 

In Eigs. 1 and 2 we show the plots of the relaxation and creep memory functions versus 
time, by assuming r = 1 and we compare them with their asymptotic representations for 
short and long times, printed with dotted and dashed lines, respectively. 
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Fig.l: The relaxation-memory function <h(t) (continuous line) with its asymptotic repre¬ 
sentations (dotted and dashed lines) versus time scaled with r. 



Fig.2: The creep-memory function 'l'(t) (continuous line) with its asymptotic representa¬ 
tions (dotted and dashed lines) versus time scaled with r. 

Finally, in Figs. 3 and 4 we show the plots of the G{t) and J{t) respectively, versus 
time, by assuming r = 1 , G'(0+) = J(0+) = 1 in Eqs. (4.18) and (4.19) and taking 100 
terms in the corresponding Dirichlet series. 



Fig.3: The normalized Relaxation Modulus G{t) versus time scaled with r. 
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Fig.4; The normalized Creep Compliance J{t) versus time scaled with r. 

Considering the asymptotic and numerical results, as well as the results in [14] and [16], 
we observe that the relaxation and creep properties of this model are consistent with those 
of a peculiar viscoelastic medium that behaves since short times as a fractional Maxwell 
solid of order 1/2 up to as a standard Maxwell solid for long times. 

6. Conclusions and final remarks 

In this paper we have shown that, from the dynamic point of view of uni-axial wave 
propagation, an elastic tube filled with a viscous fluid is equivalent to a viscoelastic solid 
with particular rheological properties that evolve as shown in Section 5. We note that 
this analogy suffers of subtleties due to the presence of Dirichlet series. In any case this 
model is worth to be investigated in future with respect to the evolution of the transient 
waves propagating through it. In fact, as formerly pointed out by [4], see also [14], the 
singular behavior of the creep function at t = 0 induces a wave-front smoothing of any 
initial discontinuity, like a diffusion effect. For more details on this smoothing effect found 
for transient waves propagating in media with a singular memory we refer to [9]. 

For that concerning the profile of the fluid inside the tube, we understand that the 
Poiseuille parabolic profile, corresponding to the Maxwell model, found for long times is 
preceded by a boundary layer profile, corresponding to the fractional Maxwell model of 
order 1/2 for short times. This fact is compatible with what expected in a pulsatile blood 
flow in large arteries for slow and high frequencies and it is essentially due to the effect 
of the blood viscosity. For further research we propose to investigate the evolution of the 
shock formation due to this boundary layer following the preliminary non-linear analysis 
by [15]. 

Finally, we note that more realistic models for ID blood flow take into account the 
viscoelastic nature of the arterial wall. On this respect the literature is wide but more 
recent articles require fractional viscoelastic models as well, see e.g. [17] and references 
therein. So we point out the necessity to include in the available ID models for the arterial 
wall viscoelasticity the effects of the blood viscosity discussed in the present paper. 
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Appendix: Mathematical discussions 
We provide the details for the proof of some statements found in the text. 


Proof of the expressions (4.3) and (4.10). Consider the Laplace Transform of the 
Relaxation-memory function 4>(s) 


( 6 . 1 ) 


5(s) = 


2 Ii(v/sr) 

a/st/o(\/st) 


where we will consider r = 1 for sake of simplicity (it can be restored by make the 
substitution s <-> sr). Then, the Relaxation-memory function is given by 


( 6 . 2 ) 


= 4 2 exp (-A^t) 


where A„ 6 M are such that Jo{Xn) = 0 , Vn = 1, 2, 3,... and t > 0. 


Proof. : Firstly, consider the power series representation for the Modified Bessel functions 
of the First Kind: 


loiVs) “ ^ + 4 + ^ + O(s^), 


one can eventually deduce that 


(6.4) 


Hs) 


2 hiVs) ^ 2 2 + B 
VsIoiVs) i + | + g + ... 


which means that the function of our concern is regular in s = 0 and it does not have 
any branch cuts. Then we can obtain the required function by means of the Bromwich 
Integral: 


(6.5) 


m 



e"* ds. 


In particular, 4>(s) has simple poles such that: Io{-\/s) = 0. Now, if we rename ^/s as 
^/s = —iX, then 

4o(V^) = 0 Jo(A) = 0. 


Moreover, 


( 6 . 6 ) 




From the previous statements, we can then conclude that: 

4>(t) = '^TZes |<I>(s)e^* 

(6.7) 


TZes 


n=l 


^ 4i(v^) St 

Vs loiVs) 


5=-A2 
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It is quite straightforward that 


oSt ' 



Io{V~s) . 

) S^Sn 

(6.8) 

= lim (s - Sn) 

S—>Sn 

2 /i(v^)e^* 


a/s Io{y/s) 


= 4 exp (snt) . 


Thus, 



(6.9) 

00 

<I>(t) = TZes |4>(s) j 

n^l 

00 

" n=l 


o-AS t 


□ 


Let us now consider the Laplace Transform of the Creep-memory function 'I'(s) 

2 hiVs) 


( 6 . 10 ) 

Then 'l'(t) is given by 

( 6 . 11 ) 


= rr ( r\- 
VshWs) 

00 

^{t) =8+ 4 2 ex.p{-filt) 

n^l 


where Hn are such that J 2 if^n) = 0 and ^ 0, for n e N. 

Proof. : By means of the same procedure shown before we are able to point out that 

^(i) = J^TZes {$(s)e"* 


( 6 . 12 ) 


TZes 


n=0 


2 hiVs) 
y/s hiVs) 


-.st 


with Hn such that J 2 (Mn) = 0, 7 ^ 0, f or n e N and hq = 0. 

Now, we have to distinguish two cases: 

If 7 ^ 0, then 

TZes 1 4 ^( 5 ) =4 exp (s„t) . 


(6.13) 

Otherwise, if Sn = po = 0, 

(6.14) 7^es|$(s)e"* 

Thus, 

00 

'l'(t) = TZes -j 'L(s) e 


= lim s 41 ( 5 ) = 8 . 

=0 s^O 


(6.15) 


n=0 


st 


4 2 exp {-fil t) 


n=l 


□ 


From the above results we are able to conclude that representation of both memory 
function <h(t) and 4'(t) are given by Dirichlet series (6.2) and (6.11), whose convergence 
is proved in the following. 
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On the convergence of the Dirichlet series (6.2) and (6.11). The series (6.2) is 
convergent for t > 0. 

Proof. : Consider a Generalized Dirichlet Series: 


(6.16) 


f{z) = 2 On exp {-OnZ) , 


Z 6 


n=l 


In general, we have that the abscissa of convergence and the abscissa of absolute conver¬ 
gence would be different, i.e. ¥= cTa, but they will satisfy the following condition: 


(6.17) 

If d = 0, then 

(6.18) 


Inn 

0 ^ iTa — CTc ^ a = lim sup- 

n—>00 On 


a = Gc = (Ta = lim sup 

n^oo 


In lunl 


O-r 


In our case On = 1 and On = ^ 0. Then, we have to understand the behavior of the 

coefficients An for n » 1, where Jo(An) = 0, Vn e N. 

Considering the asymptotic representation 


VI; 


(6.19) 

we get to the following conclusion: 

(6.20) To(An) = 0, for n » 1 
Thus, 

Inn Inn n»i Inn n^oo 


AnOCn , for n » 1. 


Olr 


( 6 . 21 ) 

which tells us that d = 0. Finally. 

( 6 . 22 ) 


AH 


n^ 


r ln|an| „ 

a = Gc = Ga = lim sup - = 0 

n-^co Oin 


being On = 1. 

This result implies that the series (6.2), with On = 1 and an = ^ *1) converges for 

Rejz} = t > 0. □ 

In a similar way a we can prove that the series (6.11) is convergent for t > 0. 

On the asymptotic representations. Finally, we derive the asymptotic representations 
for <h(t) and T(t) as t —> 0"^ applying the Tauberian theorems to the corresponding Laplace 
transforms: 


(6.23) 


$(s) = 


2 /i(V^) 


T(s) = 


2 /i(V^) 


/sr/o(V^) ’ hiV^) ' 

Then, in view of the asymptotic representation of the modihed Bessel functions as 2 ; —> oo 
with |arg(z)| < 7r/2 and for any u 

Uiz) - exp{z) , 

\ Ztt 


see e.g. [1], we conclude that for z = ^/^ cc {s > 0) we get 

$(s)~-^, T(s) ^ ^ 


ST 


ST 
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SO that 


(6.24) 




\J'KT 


^ t ^ 0+ 

V'TTT 


References 

[1] M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions, Dover, New 
York (1965). 

[2] D.W. Barclay, T.B. Moodie and V.P. Madan, Linear and non-linear pulse propagation 
in fluid-filled compliant tubes, Meccanica 16, 3-9 (2004). 

[3] A.C.L. Barnard, W.A. Hunt, W.P. Timlake, E. Varley, A theory of fluid flow in 
compliant tubes, Biophys. J. 6, 717-724 (1966). 

[4] P.W. Buchen, F. Mainardi, Asymptotic expansions for transient viscoelastic waves. 
Journal de Mecanique 14, 597-608 (1975). 

[5] N.E. Daidzic, Application of Womersley model to reconstruct pulsatile flow from 
Doppler ultrasound measurements. Journal of Fluid Engineering 136, 041102/1-15 (2014). 

[6] A. Giusti, F. Mainardi, A linear viscoelastic model for arterial pulse propagation. 
Proceedings of ICMMB-19, 168-171 (2014). 

[7] A. Giusti, F. Mainardi, On infinite series concerning zeros of Bessel functions of the 
first kind, e-print, arXiv:1601.00563 (2016). 

[8] R. Gorenflo, A. Kilbas, F. Mainardi, S. Rogosin, Mittag-Leffler Functions, Related 
Topics and Applications, Springer, New York (2014). 

[9] A. Hanyga, Wave propagation in media with singular memory, Mathematical and Com¬ 
puter Modelling 34, 1399-1421 (2001). 

[10] G.H. Hardy, M. Riesz, The General Theory of Dirichlet Series, Cambridge University 
Press, Cambridge (1915). 

[11] X. He, D.N. Ku, J.E. Moors Jr, Simple calculation of the velocity profiles for pulsatile 
flow in a blood vessel using Mathematica, Annals of Biomedical Engineering 21, 45-49 
(1993). 

[12] H.W. Hoogstraten, G.H. Smit, A mathematical theory of shock-wave formation in 
arterial blood flow, Acta Mech. 30, 145-155 (1978). 

[13] R.C. Koeller, A theory relating creep and relaxation for linear materials with memory, 
J. Appl. Mech. 77, 031008-1-031008-9 (2010). 

[14] F. Mainardi, Fractional Calculus and Waves in Linear Viscoelasticity, Imperial Col¬ 
lege Press, London (2010). 

[15] F. Mainardi, H. Buggisch, On non-linear waves in liquid-filled elastic tubes, in U. 
Nigul and J. Engelbrecht (Eds), Nonlinear Deformation Waves, Springer, Berlin, 1983, 
pp. 87-100. 

[16] F. Mainardi, G. Spada, Creep, relaxation and viscosity properties for basic frac¬ 
tional models in rheology, Eur. Phys. J. Special Topics 193, 133-160 (2011). [E-print 
http: / / arxiv. org/abs/1110.3400] 

[17] P. Perdikaris, G.E. Karniadakis, Fractional-order viscoelasticity in one-dimensional 
blood flow models, Ann. Biomed. Eng. Vol 42 No 5 1012-1023 (2014) [DOL 
10.1007/sl0439-014-0970-3] 

[18] Yu.N. Rabotnov, Elements of Hereditary Solid Mechanics, Mir Publishers, Moscow 
(1980). 




A DYNAMIC VISCOELASTIC ANALOGY 


FOR FLUID-FILLED ELASTIC TUBES 


15 


[19] S. Rogosin, F. Mainardi, George William Scott Blair - the pioneer of factional calculus 
in rheology, Commun. Appl. Ind. Math. 6, pp. 20 (2014) DOI; 10.1685/journal.cairn.481 

[20] I. N. Sneddon, On some infinite series involving the zeros of Bessel functions of the 
first kind, Proc. Glasgow Math. Assoc. 4, 144-156 (1960). 

[21] J.R. Womersley, Method for the calculation of velocity, rate of flow and viscous drag 
in arteries when the pressure gradient is known, J. Physiol. 127, 553-563 (1955). 

[22] J.R. Womersley, An Elastic Tube Theory of Pulse Transmission and Oscillatory Flow 
in Mammalian Arteries, Wright Air Development Center, Technical Report, WADC-TR, 
56-614 (1957). 

^ Department of Physics & Astronomy, University of Bologna and INFN. Via Irnerio 46, 

Bologna, ITALY. 

E-mail address: andrea.giusti@bo.infn.it 

^ Department of Physics & Astronomy, University of Bologna and INFN. Via Irnerio 46, 

Bologna, ITALY. 

E-mail address: francesco.mainardi@bo.infn.it 



